COHERENCE STABILITY AND EFFECT OF RANDOM NATURAL 
FREQUENCIES IN POPULATIONS OF COUPLED OSCILLATORS 

GIAMBATTISTA GIACOMIN, ERIC LUCON, AND CHRISTOPHE POQUET 

Abstract. We consider the (noisy) Kuramoto model, that is a population of N oscil- 
lators, or rotators, with mean-field interaction. Each oscillator has its own randomly 
chosen natural frequency (quenched disorder) and it is stirred by Brownian motion. In 
the limit N — > oo this model is accurately described by a (deterministic) Fokker-Planck 
equation. We study this equation and obtain quantitatively sharp results in the limit 
of weak disorder. We show that, in general, even when the natural frequencies have 
zero mean the oscillators synchronize (for sufficiently strong interaction) around a com- 
mon rotating phase, whose frequency is sharply estimated. We also establish the stability 
properties of these solutions (in fact, limit cycles). These results are obtained by identify- 
ing the stable hyperbolic manifold of stationary solutions of an associated non disordered 
model and by exploiting the robustness of hyperbolic structures under suitable pertur- 
bations. When the disorder distribution is symmetric the speed vanishes and there is a 
one parameter family of stationary solutions, as pointed out by H. Sakaguchi [20] : in this 
case we provide more precise stability estimates. The methods we use apply beyond the 
Kuramoto model and we develop here the case of active rotator models, that is the case 
in which the dynamics of each rotator in absence of interaction and noise is not simply 
a rotation. 
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1. Introduction 

1.1. Collective phenomena in noisy coupled oscillators. Coupled oscillator models 
are omnipresent in the scientific literature because the emergence of coherent behavior 
in large families of interacting units that have a periodic behavior, that we generically 
call oscillators, is an extremely common phenomenon (crickets chirping, fireflies flashing, 
planets orbiting, neurons firing,...). It is impossible to properly account for the literature 
and the various models proposed for this kind of phenomena, but while a precise descrip- 
tion of each of the different instances in which synchronization emerges demands specific, 
possibly very complex, models, the Kuramoto model has emerged as capturing some of 
the fundamental aspects of synchronization [I]. It can be introduced via the system of N 
stochastic differential equations 

K N 

d<fi(t) =Ujdt--Y, " d* + adB j (t) , (1.1) 

i=i 



for j = 1, . . . , N , where 
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(1) {-Sj},- = i is a family of standard independent Brownian motions: in physical 
terms, this is a thermal noise; 

(2) {wj}j=i—N is a family of independent identically distributed random variables of 
law [i: they are are the natural frequencies of the oscillators and, in physical terms, 
they can be viewed quenched disorder; 

(3) K and a are non-negative parameters, but one should think of them as positive 
parameters since the cases in which they vanish have only a marginal role in the 
what follows. 

The variables (p" are meant to be angles (describing the position of rotators on the circle 
S), so we focus on ipj mod 2ir and (jl.ip defines, once an initial condition is supplied, a 
diffusion process on S . Note that if {^(-)}j=i,...,iV solves (jl.ip . also {iffj (•) + <p}j=i t ... t N , 
with (/? G S, is a solution: this is the rotation symmetry of the system that will repeatedly 
make surface in the remainder of the paper. 

Some of the main features (jl.ip are easily grasped: each oscillator rotates at its own 
speed, it is perturbed by independent noise and it interacts with all the other oscillators: 
the interaction tends to align the rotators. It may be helpful at this stage to point out 
that if [i = Sq, that is the natural frequencies are just zero, then the dynamics is reversible 
with invariant probability measure that, up to normalization, is 



where A at is the uniform measure on §> N . The Gibbs measure in (|1.2p is a well known 
statistical mechanics model - it is the classical XY spin mean field model or rotator mean 
field model - treated analytically in |19j [T7] in the N — > oo limit. In particular, the 
model exhibits a phase transition at K = K c := l/a 2 , that is effectively a synchronization 
transition: in the N — > oo limit we have that for K < K c the rotators become independent 
and uniformly distributed over S, while for K > K c the limit measure is obtained by 
choosing a phase 9 uniformly in § and by choosing the values of the phase of each oscillator 
by drawing it at random following a suitable distribution that concentrates around 9. 
However, in [3J Prop. 1.2], it is shown that, unless fi = 5q, the model is not reversible 
(for \x almost surely all the realization of uj) and one effectively steps into the domain of 
non-equilibrium statistical mechanics. 

Our approach actually relies on a sharp control of the reversible case and works when 
the system is not too far from reversibility, that is for weak disorder. Our approach 
actually applies well beyond (II. ID : here we will treat explicitly the case ojj is replaced by 
U(tffj,ojj), that is the natural frequency Uj is replaced by a natural dynamics that can be 
substantially different from one oscillator to another. This model is a disordered version 
of the active rotator model considered for example in [21] . 

Since we will focus on a > 0, from now on, for ease of exposition, we set a := 1. 

1.2. The Fokker-Planck or McKean-Vlasov limit. An efficient way to tackle (jl.ip 
is to consider the empirical probability on S x R 




(1.2) 




(1.3) 
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In fact, in the N — > oo limit, the sequence of measures {vfj t } jv=l,2,... converges to a 
limit measure whose density (with respect to Ai ® fi) solves the nonlinear Fokker-Planck 
equation 

d tPt (e,u) = ^Ap t (e,u)-de( Pt (e,u J )((j*p t )^6) + u J )y (i.4) 

where J(8) = —K sin(#), * denotes the convolution and (•}„ is a notation for the integration 
with respect to /i, so (J * u)^(9) = f R L J{tp)u{9 — 93, u) dtpfi( du) is the convolution of J 
and u, averaged with respect to the disorder. Here and throughout the whole paper A 
means d%. The Fokker-Planck PDE (jl.4p appears repeatedly in the physics and biology 
literature, see e.g. [TJ [212 122], and a mathematical proof (and precise statement) of the 
result we just stated can be found in [5j[l3]. Notably, in [13] the result is established under 
the assumption that J \u\fi( doj) < 00 and emphasis is put on the fact that the result holds 
for almost every realization of the disorder sequence {wjj^i^,...- Let us point out that in 
(jl.4p a; is a one dimensional real variable, while in (jl.ip the superscript w is a short for the 
whole sequence of natural frequencies. Since what follows is really about (11. 4p this abuse 
of notation will be of limited impact. 

In AppendixEl we detail the fact that (ll.4p generates an evolution semigroup in suitable 
spaces. Here we want to stress that (ll.4j) can be viewed as a family of coupled PDEs, one for 
each value of us in the support of fi: pt(-,ui) is the distribution of phases in the population 
of oscillators with natural frequency uj. 

1.3. About stationary solutions to (jl.4p . Remarkably ([20], see also [H]), if /i is 
symmetric all the stationary solutions to (|1.4p can be written in a semi-explicit way as 
q(9 + 9q,u) (9q is an arbitrary constant that reflects the rotation symmetry) where 

S{9,u,2Kr) 

q V' u) := Z(u,2Kr) ' (L5) 



with 



S(0,u,x) = e^ 9 '^ 



[1 - e 4 ™) / e~ G( - u '^ du + e 4 ™ / e 
Jo Jo 



(1.6) 



and G(u,y,x) = xcos(u) + 2yu, Z(co,x) = J s S(9,uj,x)d9 is the normalization constant 
and r S [0, 1] satisfies the fixed-point relation 

„, s f Lcos(9)S(9,uj,x)d9 

r = V(2Kr), where V(x) := , V — fi(du) . (1.7) 

J R Z(u,x) 

A series of remarks are in order: 

(1) r = solves (|1.7p and this corresponds to the fact that q(-) = ^7 is a stationary 
solution. It is the only one as long as K does not exceed critical value K c which 
is in any case not larger than 



._ l + 4w 2 , 

as one can easily see by computing (see e.g. |llj ) the derivative of ^^(2^-) at 
the origin and noticing that is larger than one if and only if K > K and that 
'f(-) < 1, see Figured) 
(2) When (jl.7j) admits a fixed point r > 0, and this is certainly the case if K > K, a 
nontrivial stationary solution is present and in fact, by rotation symmetry, a circle 
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(A) Function ^»(2K-) for (b) Function ^(2K-) for 

V = 5 ,K = 2. y, = |(tf_i + «i), if = 

3.5. 



Figure 1. Plot of the fixed-point function * M (2A'-) for two choices of K and p. * 4 °(-) is 
strictly concave with derivative at the origin equal to 1 /2 (Fig. I lap but even for a simple 
instance of fj, (Fig. Ilb|) concavity is lost and there are several non-trivial fixed-points, 
each of them corresponding to one circle of non trivial stationary solutions. Note that 
in the case of Fig. lib! K < K = 5 so that the phase transition is not given by the 
derivative of ty^(2K-) at the origin. 

of non-trivial stationary solutions. Such solutions correspond to a synchronization 
phenomenum, since the distribution of the phases is no longer trivial. 

(3) As explained in Figure Q] and its caption, in general there can be more than one 
fixed point r > 0: in absence of disorder there is only one positive fixed point 
(when it exists, that is for K > 1), but this fact is non-trivial even in this case (see 
below). Uniqueness is expected for fx which is unimodal, but this has not been 
established. 

(4) While the local stability of is understood [22] and it holds only if K < K, the 
stability properties of the non-trivial solutions are a more delicate issue. 



1.4. An overview of the results we present. Here are two natural questions: 

• What are the stability properties of the non-trivial stationary solutions? 

• What happens if \i is not symmetric? 

Our work addresses these two questions and provides complete answers for weak disor- 
der. The precise set-up of our work is better understood if we remark from now that we 
can assume := j dw) = 0. In fact, if this is not the case we can map the model to a 
model with m u = by putting ourselves on the frame that rotates with speed m u , that is 
if we consider the diffusion {^(t) — m^i}^!,...^- So, we assume henceforth ra w = and 
we rewrite the natural frequencies as 5uj, with 5 a non-negative parameter. We assume 
moreover that 



Supp(^) C [-1,1] 



(1.9) 
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In this set-up, (|1.4|) becomes 

d tP s t (9,u) = ±Ap*(e,u;)-d e (p s t (0,oj)((J*pi}^e) + 5uj)y (1.10) 

Note that this leads to (obvious) changes to fjl .5j) - f)l . 7[) . We have introduced this param- 
eterization because the results that we present are for small values of S. In particular we 
are going to show that for any K > 1, there exists 5q > such that for S G [0, 5q] 

• there exists a solution p^(9,uj) to f 1 1 . 1 j) of the form q(9 — c fl (5)t), we show that 
0^(5) = 0(5 3 ) and we actually give an expression for \iras\o c fM (5) / 5 3 : this is a 
rotating wave (or limit-cycle) for the dynamical system (jl.lOp and we establish its 
stability under perturbations; 

• when /i is symmetric and K > K we show that there is, up to rotation symmetry, 
only one non-trivial solution and that it is (linearly and non-linearly) stable. 

The results we obtain are based on the rather good understanding that we have of the 
case (5 = that, as we have already explained, is reversible and the corresponding Fokker- 
Planck PDE is of gradient flow type (e.g. [15] and references therein). These properties 
have been exploited in [3] in order to extract a number of properties of the Fokker-Planck 
PDE (denoted from now on: reversible PDE) 

d tPt {Q) = l -Ap t {9)-de{pt{9){J*ptW)), (1.11) 

and notably the linear stability of the non-trivial stationary solutions. In fact one can find 
in [3] an analysis of the evolution operator linearized around the non-trivial stationary 
solutions. Some of the results in [3] are recalled in the next section, but they are not 
directly applicable because the 6 = case that corresponds to what interests us is rather 

d t p t (9,co) = ±Apt(0,u>) -d e (p t (9,u)((J*p t }^9^ , (1.12) 

which we call non-disordered PDE. So the natural frequencies have no effective role beyond 
separating the various rotators into populations with given natural (ineffective) frequency 
that now are just labels. But in order to set-up a proper perturbation procedure we need 
to control (| 1 . 12[) and, in particular, we need (and establish) a spectral gap inequality for 
the evolution (|l,12p linearized around the non-trivial solutions. 

This spectral analysis is going to be central both for the general and for the symmetric 
disorder case. In the general set-up we are going to exploit the normally hyperbolic 
structure [10\ [18] of the manifold of stationary solutions of (|1.12p and the robustness of 
such structures (like in [8] ) . In the case of symmetric \x we can get more precise results by 
ad hoc estimates, made possible by the explicit expressions (|1.5p - (jl.7p . and use results in 
the general theory of operators |16] and perturbation theory of self-adjoint operators |12] , 

The normal hyperbolic manifold approach allows to treat cases that are substantially 
more general and notably the case of 

d tPt (9,u) = ^A Pt (9,u;) - d e (pt(9,u)((J * p t )v(9) + 5U(9,u))y (1.13) 

which is the large N limit of (jl.ip with the term ujj dt replaced by U(ip l j (t),ujj) dt, with 
U G C^SxI; R). In this case each oscillator has its own non-trivial dynamics which may 
be very different from the dynamics of other oscillators: consider for example 



U(cp,U}) = b + oj + a sin((/?) , a, b € R , 



(1.14) 
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and fi uniform over [—1, 1]. For a £ (—1, 1) there are some active rotators |21|, [3] that in 
absence of noise and interaction (cr = K = 0) rotate (this happens if \b + u\ > \a\ and of 
course the direction of rotation depends on the sign of b + w) and others that instead are 
stuck at a fixed point (this happens if \b + w| < \a\). Our approach allows us to establish 
that there is a synchronization regime for K > 1 and 5 small and to describe the dynamics 
of the system in this regime. This is going to be detailed in Section [5j 

The two questions raised at the beginning of this section have been already repeatedly 
approached but looking at synchronized solutions as bifurcation from incoherence. The 
results are hence for K close to the critical value corresponding to the breakdown of linear 
stability of \/2tt: one can find a detailed review of the vast literature on this issue in [H 
Sec. III]. Our results are instead for arbitrary K > 1, but 5 smaller than 5q{K) and of 
course 5q(K) vanishes as K approaches 1. 

2. Mathematical set-up and main results 

2.1. The reversible and the non-disordered PDE. We first recall some results about 
the reversible PDE (jl.lip . The stationary solutions qo(6) = q(9,0) are, up to rotation 
invariance, given by (|1.5p - (|1.7p . but formulas get simpler, namely 

90 W = 7 (OK \ eX P( 2 ^ COS(0)) , (2.1) 

where Zo(x) := Z(0, x)a and this time we have the more explicit expression Zq(x) = 
(x) is the normalization constant and ro is a solution of the fixed- 
point problem 

r = y {2Kr ) where # (x) := yj^ , (2-2) 
where we used standard notations for the modified Bessel functions 

Ii(x) = — [(cos(6)yexp(xcos(9))d6 i = 0, 1 . (2.3) 

The mapping is increasing, concave (see p2]) and with derivative at equal to i. 
Consequently if K ^ 1, ro = is the unique solution of the fixed-point problem, and 
q(-) = j- is the only stationary solution of (jl.lip . If K > 1, we get in addition a circle 
(because of the rotation invariance) of nontrivial stationary solutions 

M rev := {g^o(-) := g (- - V) : V G §} with q (9) := ex P( 2 ^o cos(g)) 

J s exp(2A ro cos(0)J 



where ro = ro(K) is the unique non trivial fixed-point (12. 2p . 

Let us now focus on the non-disordered PDE (| 1 . 12[) and let us insist on the fact that we 
are interested in solutions such that oj ) is a probability density. Observe then that if 
q(6,uj) is a stationary solution of (|1.12p . we see (Appendix |A]) that q is C°° with respect 
to 9 and that (q)^ is a stationary solution for (jl.lip . So there exists ip G S such that 
= Qtp an d a short computation leads to 

(J*q)^9) = -Ksm(9-^P), (2.5) 

and, since Jgq(6,uj)d9 = 1 for almost all w, we obtain that q(-,uj) = q^j(-) for almost 
all u. In conclusion, with some abuse of notation, we can say the stationary solutions of 
(jl.lip and (jl,12p are the same: of course in the second case the function space includes 
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I i|2 II 1 1 CO || ||2 (1 Cl\ 

Fll-i.fei,^ ^ TfTli H tt H-i,fc2,At ' ^ y ) 



the dependence on u, so we choose a different notation, that is Mq, for the corresponding 
circle of non-trivial stationary solutions. 

An important issue for us is the stability of Mq (for its existence we are assuming K > 1) 
and for this we denote by A the linearized evolution operator of (|1 . 12[) around go 

Au(9,u) := ±Au(e,u})-d 9 (qa(e){J*u)r(e)+u(0,u)J*qo(oj) (2.6) 
with domain 

T>(A) := lu G C 2 (§ xR,R) : J u(9,u) d0 = for all u\ . (2.7) 

For any smooth positive function k : 8 i— > R, we introduce the Hilbert space defined 
by the closure of T>(A) for the norm || • ||-i,fc )jU associated with the scalar product 

where w a.s., U(-,ui) is the primitive of u(-,w) such that J* § ^jjjy- d# = 0> an d V(-,u;) is 
defined in the analogous fashion. Let us remark (see [SJ Sec. 2]) immediately that 

IMoo || 1 1 2 

N 

so that all the norms we have introduced are equivalent. For the case k{ ) = 1 we use the 
notations H' 1 and || • H-i,^. We will prove the following result, which is just technical, 
but it will be of help to understand our main results: 

Proposition 2.1. A is essentially self-adjoint in H~ „. Moreover the spectrum lies in 
(—oo,0], is a simple eigenvalue, with eigenspace spanned by dgqo, and there is a spectral 
gap, that is the distance Xk between and the rest of the spectrum is positive. 

The proof of this result builds on [3j Th. 1.8] that deals with the reversible case and the 
(lower) bound on the spectral gap Xk that we obtain coincides with the quantity X(K) 
in [31 Th. 1.8] (this bound can be improved as explained in in [31 Sec. 2.5] and sharp 
estimates on the spectral gap can be obtained in the limit K \, 1 and K /*■ oo). For the 
reversible evolution, the linear operator L qo is defined by 

L qo u{0) ■= ~ de (qo(O)J * u{9) + u{6)J * q (9)^ , (2.10) 

with domain D{L qo ) given by the C 2 (S,1R) functions with zero integral. 

2.2. Synchronization: the main result without symmetry assumption. Proposi- 
tion [2J] is a key ingredient for our main results and the functional space H^ 1 appears in 
it, but an important role is played also by L 2 (A ® fx), X is the Haar measure on S, whose 
norm is denoted by || • ||2 lAt . For C > and M C L 2 (A®/x) we set A/2 lAt (M, C) := {u : there 
exists v G M such that \\u — f H2 u < C}. In the statement below q G Mq is the element of 
the manifold such that q(-,uj) = qo(-), cf. (|2.ip . with ro(K) > (hence K > 1). 

Theorem 2.2. For every K > 1 there exists 5q = 5q{K) > such that for \8\ < 6q there 
exists qs G L 2 (X ® fj), satisfying \\qs — q\\2,^ = 0(5) and a value c M (<5) G M such that if we 
set 

q?\9,u) := &(0-c„(5)t-^), (2.11) 
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then q^ solves (jl.lOp . Moreover 

(1) the family of solutions {q^}^ is stable in the sense that there exist two posi- 
tive constants j3 = f3(K) and C = C{K) such that if Pq £ A/2,u(-Mqj $)> an d 
J s Pq(8,lo) d9 = n(duj)-a.s., then there exists £ § such that for all t > 

Wq^-pfhn < 2Cexp(-/3t). (2.12) 

(2) we have 

c,(5) = ^ 9 rl + 0(5 5 ) , (2.13) 

where is the unique solution of 

An {2) = cod e n W and (n {2 \d e q ) = 0, (2.14) 

and nfl' is the unique solution of 

An W = uj dgq and (n (1) , d e q Q ) =0. (2.15) 

In the proof of Theorem 12.21 one finds also further estimates, in particular (see (|4.19p ) 
that one has 

q s = q + <5n (1) + 5 2 n {2) + L 2 (5 3 ) . (2.16) 
Actually, see Remark l4.2l the argument of proof can be pushed farther to obtain arbitrarily 
many terms in development (|2,16p . as well as in 

c»($) = c 3 6 3 + c 5 5 5 + . . . . (2.17) 

In Table [T] we report a comparison between the c„(5) obtained by solving numerically 
(jl.lOp and by evaluating the leading order i.e. by using f|2. 13jl . 



s 


K = 2 


K = 1.5 


K = 1.1 


0.5 


-1.56300 • 10" 2 


-8.59626 • 10~ 2 


-3.01064 • 10" 1 


0.1 


-1.23998 • 10~ 2 


-6.84835 • 10~ 2 


-2.72117 • 10" 1 


0.05 


-1.23072 • 10~ 2 


-6.79553 • 10~ 2 


-2.69460 • 10" 1 


0.01 


-1.22776 • 10~ 2 


-6.77921 • 10~ 2 


-2.68603 • 10" 1 


0.005 


-1.22767 • 10~ 2 


-6.77869 • 10~ 2 


-2.68576 • 1Q- 1 



C3 


-1.22764 • 10~ 2 


-6.77851 • 10~ 2 


-2.68567 • 10- 1 



Table 1. For the case /i = p8i- p + (l—p)5- p , p — 0.2, we have computed (numerically) 
c M (J)/5 3 for three values of K and five values of 8. In the last line we report the value 
C3 = lim^o c ^(8)/8 that one obtains by using (|2.13[) . 



2.3. Symmetric disorder case. Let us focus on the case in which the distribution of 
the disorder fx is symmetric. In this case, at least for small disorder, Theorem 12.21 is just 
telling us that the leading order in the development for the speed Cu(S) is zero: one can 
actually work harder and show that such a development yields zero terms to all orders. In 
reality in this case we already know, see (|1.5p - (|1.7p . that for K sufficiently large there is 
at least a non-trivial stationary profile, hence, by rotation symmetry, at least one whole 
circle of stationary solutions. Actually, we can show that for 6 small there is just one 
circle, that we call Mg, of non-trivial stationary solutions and this circle converges to Mq 
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as 8 \, (in C 3 , for every j) so the rotating solutions found in Theorem 12.21 must be the 
stationary solutions in Ms- 

In order to be precise about this issue, we point out that (jl.5p - (|1.7|) are written for (|1.4|) 
while we work rather with (ll.lOp . The changes are obvious, but we introduce a notation 
for the analog of (|1.7|) : 

, a, , f L cos(6) S (9 , 5uj , x) dd 
r s = ^(2Kr s ), where, ^(x) := / ™ K ' K ' p(du) . (2.18) 



Lemma 2.3. For all K m \ n < K max , there exists 8\ = 8i(K m i n , K ma _ x ) > such that, for 
all < K m { n < K < K max and all 8 ^ 5\ the function ^ is strictly concave on [0, 1]. 
Therefore for (|1.7|) has only a positive solution rg = r§(K,p). Moreover lim^o 7 ^ = tq. 



We point out that in spite of the fact that ^ is explicit (cf. (|2.2p ). it is not so 
straightforward to show that it is concave. We show that tyg remains strictly concave for 
a small 8 via a perturbation argument. But the conjecture (see |11| and [5]) that VE^ is 
strictly concave for unimodal distributions p is still an open issue. 

Remark 2.4. A direct computation shows that the derivative of ^ at the origin is 
1/(2K 5 ), for K S := (f R 

i+4pl 2 ) course K\ coincides with K, introduced in (jL£J)). 

Under the hypothesis of Lemma 12.31 one therefore sees that there is a synchronization 
transition at K = K$ in the sense that for K < K$ the only stationary solution is ^~ while 

for K > K$ also the manifold of non-trivial stationary solutions appears (and there is no 
other stationary solution). 

Theorem 12.21 provides a stability statement for Mg. This result can be sharpened and 
for this let us introduce the linear operator 

L%u{0,u) := ^Au(e,u})-d g (u(e,u)({J*q)r(e) + 6u) + q(6,6w){J*u) li (0)), (2.19) 

The domain T>(L^) of the operator is chosen to be the same as for A, cf. (|2.7p . 

We place ourselves within the framework of Lemma 12.31 m the sense that 5 is small 
enough to ensure the uniqueness of a non-trivial stationary solution (of course existence 
requires K > K§ and this is implied by K > 1 if 8 is sufficiently small). We prove a 
number of properties of the linear operator (|2.19p . saying notably that it has a simple 
eigenvalue at zero and the rest of spectrum is at a positive distance from zero and it is in a 
cone in that lies in the negative complex half plane. We summarize in the next statement 
the qualitative features of our results on , but what we really prove are quantitative 
explicit estimates: the interested reader finds them in Section [6) 



Theorem 2.5. The operator has the following spectral properties: is a simple eigen- 
value for Lg, with eigenspace spanned by (9,u) h-» q'(9,uj). Moreover, for all K > 1, 
p G (0,1), a G (0, 7r/2), there exists 82 = 82(K, p,a) such that for all ^ 8 ^ 82, the 
following is true: 

• Lg is closable and its closure has the same domain as the domain of the self-adjoint 
extension of A; 
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• The spectrum of lies in a cone C a with vertex and angle a 

C a := |a G C ; | + a ^ arg(A) < y - a j C {z G C ; ^ 0} ; (2.20) 

• There exists a' G (0, ^) suc/i that is the infinitesimal generator of an analytic 
semi- group defined on a sector {A G C, | arg(A)| < a'}; 

• The distance between and the rest of the spectrum is strictly positive and is at 
least equal to p\x, where Xk is the spectral gap of the operator A introduced in 
Proposition \2.1\ 

2.4. Organization of remainder of the paper. In Section [3] we introduce the notion 
of stable normally hyperbolic manifold, we recall its robustness properties, and show that 
Mq is in this class of manifolds. The essential ingredient is Proposition 12. II that, directly or 
indirectly, plays a role in each subsequent section. Section [3] is also devoted to the proof of 
Proposition 12.11 The proof of Theorem 12.21 is then completed in Section U that is mainly 
devoted to perturbation arguments. The case of the active rotators is treated in Section [5l 
while Section [6] deals with the case symmetric disorder distribution and, notably, with the 
proof of Theorem 12.51 and of a number of related quantitative estimates. 



3. Hyperbolic structures and periodic solutions 

In this section we present the arguments proving the existence of the periodic solution 
of Theorem 12.21 We rely on the fact that the circle of stationary solutions Mq is a stable 
normally hyperbolic manifold, and on the robustness of this kind of structure : adding the 
perturbation term —5d$(pt(9,uj)uj) in (j!.12j) . this manifold Mq is deformed into another 
manifold M$, and thanks to the rotation invariance of the problem, Ms is a circle too. 
The spectral gap of operator A (Property 12. ip which induces the hyperbolic property of 
Mq is proved at the end of this section. 

3.1. Stable normally hyperbolic manifolds. We start by quickly reviewing the notion 
of of stable normally hyperbolic manifold (SNHM). The evolution of (jl.lOp will be studied 
in the space defined by 

Xl := ju G L 2 (\®fi), J u{e,u)&Q = 1 wa.s.j (3.1) 

where A denotes the Lebesgue measure on §. This is made possible by the conservative 
character of the dynamics. The L 2 -norm with respect to the measure A<8>// will be denoted 
by || • 1 1 2,/j, - We will also use the space X® defined by 

X° := G L 2 (\®n), J u(0,u))d6 = u a.s.j . (3.2) 

To define a SNHM, we need a dynamics: we have in mind (jl.lOp but for the mo- 
ment let us just think of an evolution semigroup in X^ that gives rise to {ut}t>o, with 
uq = u, to which we can associate a linear evolution semigroup {&(u, t)} t >o in satis- 
fying dt&(u,t)v = L(t)<&(u,t)v and &(u,0)v = v, where L(t) is the operator obtained by 
linearizing the evolution around Ut- 

For us a SNHM M C X^ (in reality we are interested only in 1-dimensional manifolds, 
that is curves, but at this stage this does not really play a role) of characteristics Ai, A2 
(0 < Ai < A2) and C > is a C 1 compact connected manifold which is invariant under 
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the dynamics and for every u G M there exists a projection P°(u) on the tangent space 
of M at u, that is lZ(P°(u)) =: T U M, which, for v G Lq, satisfies the following properties: 

(1) for every t > we have 

${u,t)P°(u )v = P°(u t )<S>(u,t)v, (3.3) 

(2) we have 

||$(u,t)P°(uoM2,„ < Cexp(A 1 t)|| t ;|| 2 , At , (3.4) 
and, for P s := 1 — P°, we have 

\mu,t)P s {u )v\\ 2 ^ < Cexp(-X 2 t)\\v\\ 2 , fl , (3.5) 

for every t > 0; 

(3) there exists a negative continuation of the dynamics {u t }t<o and of the linearized 
semigroup {&(u, t)P°(uo)v}t<o and for any such continuation we have 

||$(n,t)P o (n ) V || 2 ^ < Cexp(-Ait)|M| 2j/i , (3.6) 

for f < 0. 



3.2. Mo is a SNHM. First of all: the dynamics on Mq is trivial. For q^ G Mo, the 
projection P°^ on the tangent space is the projection on the subspace spanned by q'^: 



p ^ u = . . _i^ g ; (3.7) 



and since the dynamic on the manifold is trivial, we are allowed to choose for the param- 
eters Ai = and A 2 = Xk (where we recall that Xk is given by Proposition 12. ip . 

We are in the same situation as in [8]. For a suitable perturbation and if 5 is small 
enough, the circle Mo is smoothly transformed into another SNHM M§, which is close to 
Mq. The proof is the same as in [HI Sec. 5], which, in turn builds on results in [IS]): the 
spaces we are working in are more general since we have to deal with the disorder. Here 
suitable perturbation means being an element of C 1 (X'j > l , H^ 1 ), but it is clearly the case 
for the perturbation u t— >■ —5ujdgu when /j, is of compact support. The following theorem 
works for all C l {X^,H~ l ) perturbations: 

Theorem 3.1. [SJ Sec. 5] For every K > 1 there exists do > such that if 5 £ [0, 5q] there 
exists a stable normally hyperbolic manifold Ms in X^ for the perturbed equation (11.101) . 
Moreover we can write 

M 5 = {q^ + <f>s M (3.8) 
for a suitable function <ps G C 1 (Mo,X°) with the properties that 

• there exists C > such that sup^d^ (q^) || 2jAt + {{d^siq^h^) < C5. 



Remark 3.2. A byproduct of the proof in [SJ Sec. 5] is also that Ms is the unique 
invariant manifold in a L 2 (X, /^-neighborhood of Mo. So in the case of (jl.lOp . thanks 
to the symmetry of the problem that tells us that any rotation of Ms is still a invariant 
manifold, Ms is in fact a circle, and that the dynamics on this circle is a traveling wave of 
constant (possibly zero) speed c^(5). So the invariant manifold we get for (jl.lOp is even 
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C°°. In this sense, when dealing with (|1.1U|) . we are using only part of the strength of 
Theorem 13. II Of course this symmetry argument does not apply when dealing with (jl.13)) . 



Remark 3.3. Theorem 13.11 addresses the existence and the linear stability of the manifold 
Mg. The non-linear stability statement in Theorem 12.2( 1) follows from Theorem 13.11 
combined with [9j Theorem 8.1.1], when the dynamics is periodic with non zero speed on 
Ms- If Mg is a manifold of stationary points, the argument for the non-linear stability 
follows by repeating the argument in [TJ Th. 4.8], where the non-disordered case is treated. 

We now prove Proposition 12.11 and thus that Mq is a SNHM. 

3.3. The spectral gap estimate (proof of Proposition [27T]) . We start by remarking 
that A is symmetric for the scalar product (• , qo (recall (|2.8p ). In fact, for u and v in 
T>(A), a short computation gives (in the following we use the notation u'(9, uj) = deu(9, uj)) 

'V(0,u) fu'(9,uj) 



(v, An)_ 1 



xS 



qo (9) V 2 
1 f u(9,lj)v(9,oj) 



u(9,aj)J*qo(9)-q (9){J*u) ll (9) 



d9dn 
(3.9) 



, d9d/i + v(6,u)J *u(9,uj')d9dfi(g) fi, 

2 JRxS Qo{9) JrJ(S) 2 

where J (9) = Kcos(9). We now first prove an inequality for A that is stronger than the 
spectral gap inequality and then deduce that A is (essentially) self-adjoint. We define the 
two following scalar products, which were used for the non-disordered case in [3]: 

where U{-) is the primitive of u(-) such that J* g d9 = and 

We denote the closures of V(Lq) for these scalar products respectively by H q ^ and L qQ . 
In the disordered case, L 2 qo corresponds to the space L;L „, which we define by the closure 
of T>(A) with respect to the norm || • ||2,<j ,m associated with the scalar product 

f f u(9,uj)v(9,uj) , . 

Jr Jn qo{U) 

The two Dirichlet forms for the disordered and non-disordered case are respectively 

= ~{Au, u)_ lm ^ , (3.13) 

and 

£(u) = -(L qo u,u)_ lqo . (3.14) 

As in [3], we first prove a spectral gap type inequality that involves the scalar product 
(■ , -)2 qo - For this we introduce the projections on the line spanned by q' in the spaces 

L 9„,M and L lo 

P 2 , q0tfi u = ^] %) ^ q' for aU« = u(9,u) G , (3.15) 

WO ) %/2,q 
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and ^ 

P2, qo u = l^^J^ o for all uGLl. (3.16) 

Remark that since q' does not depend on u, 

Wo . 9o>2,«,,/i = Wo > 9o> 2 , go and (q' , g&)_ 1>W)>/i = (<?(, , <Zo)_i, go , (3-17) 
and that for all u S i/!„ ,, 

ft,«o,M u = ( P 2,qo u )h = p 2, go (u)n- (3-18) 
Proposition 3.4. For all u £ L qo ^ such that for almost every uj, J^u(-,lj) = 

S^iu) > c A - (« - i^ ia o )A4 « , u - P2, qo ^u) 2 qo li , (3.19) 

with 

c K = 1 - K(l - rl) e (0, 1/2) . (3.20) 

The proof of this proposition relies on the corresponding result for the non-disordered 
case.: 

Proposition 3.5. (see [3, Prop. 2.3],) For all u G Lg suc/i that for almost every u, 
f s u(-,u) = 

£ (v) ^ c K (u- P 2 , qo u , u - P2, qo u) 2m • (3.21) 

Proof of Proposition (|3.4p . The first step of the proof is to make the Dirichlet form of the 
non-disordered case appear in the the disordered case one, that is 

S^iu) = (£(u)} M + / / u(O,u)J*[u(0,u)-u(6,u/)]d6dn®ii, (3.22) 
Jm. J(s) 2 

= (£(u)) f , + - [ [ [u(9,u)-u(6,uj')}J*[u(9,uj)-u(9,uj')}d8dfi®fi, (3.23) 
2 Jr J(S) 2 

and from Proposition (|3.5p we see that 

(£ (u))^ > C K (u - P2, qo U , U - P2,q u) 2 qo . (3.24) 

Now remark that if we define 

V = U - P2,q ,fiU , (3.25) 

using (|3.18p we get 

V - P2,q V = U - P2,q U , (3.26) 

and so 

{£(u))p ^ c K (v- P 2 , qo v, v - P2, qo v) 2 qo41 • (3.27) 

We now introduce an orthogonal decomposition of the space L? which is well adapted 
to the convolution with J. 

Lemma 3.6. (See [31 Lemma 2.1]. ) We have the following decomposition 

L l = F o ® X F l/2 ® X Fk-i/2 (3-28) 

where 

^ ao + ^2 a i cos 0^) + b i s[n (j8) ; °i + < 00 ( ( 3-29 ) 

3 > 2 j I 
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and both F 1 / 2 and F K _i/ 2 are one dimensional subspaces generated respectively by ^ 
s'm(8)q(8) (= —q' Q (9) /2Kro) and by 9 ^ cos(9)q (9). Moreover, when u G F\, then 

J* u = —u. (3.30) 
Qo 

With the help of Lemma 13.61 we can find a lower bound for the last term in (|3.23p : 
choose a such that P2, qo u = aq' , so that we can write 



Wo > q'o)2, qo 



£n(u) ^ c K (v-P 2 , qo v, v- P 2 , qo v) 2 + — '— / (a(w)-a(w)) d/x®/x. (3.31) 

But if Pz,q v = (3q' (recall that v = u — P2,q ,^u), then since P2, qo , l iU is colinear to q' , for 
almost all uj, 



P(u) - P(<J) = a(uj) - a(oj') (3.32) 
and since v is orthogonal to q' (with respect to (•, -)2, qo ,^) we get 

/3(w)d// = 0. (3.33) 

So (I3.3ip becomes 

W > C K (v ~ P 2 , qo v , v - P2, qo v) 2 ^ + i% ' f 2 '" J lV) d^ • (3.34) 

It is sufficient to compare this last minoration with the norm (v , v) 2 qo , and from Lemma 
13.61 it comes 

( v > U ) 2 ,< ?0 ,M = ( v - p 2,<ioV, v ~ P2, qo v) 2 ^ + (q' , q' ) 2 qQ J P 2 (uj)dfi. (3.35) 

This completes the proof of Proposition 13.41 □ 
We now need two lemmas comparing the scalar products (• , -) 2 9o „ and (• , ■)_i gQ „. 
They correspond to Lemmas 2.4 and 2.5 in [3]. Their proofs are very similar to the 
proofs of the results corresponding results in [3] (to which we refer also for the explicit 
values of the constants C and c appearing below) and they use in particular the rigged 
Hilbert space representation of H~ „ (see [U p.82]): namely, one can identify Hz^ as 
the dual space V of the space V closure of V(A) with respect to the norm := 



Jmx§ v '{®i UJ ) 2 d^A*( dw) J . The pivot space H is the usual L 2 (X (g) ^) (endowed with the 



Hilbert norm ||u||2,/i := ( Jr x § w) 2 d#//( dcj)J ). In particular, one easily sees that 

the inclusion V C H is dense. Consequently, one can define T : H — > V' by setting 
Tu(v) = J RxS u(9,uj)v(9,uj) d9n(dco). One can prove that T continuously injects -ff into 

V' and that T{H) is dense into V' so that one can identify u G H with Tu G V'. Then 
for u G H, 

rp II JW rju? 

\\u \\v> = U u \\v = sup -77— 7. — = \ / — , (3.3b) 
v£V \\v\\v V J q ° 

which enables us to identify H~ Q K with V' . 
We define the projection in H~ Q •. 

( u , %) -l,q ,n i fQ „„x 

P-i, qo ,^u = £q . (3.37) 

WO ' %/-l, qo 
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Lemma 3.7. For every K > 1 there exists a constant C = C(K) > such that for u G L 2 
such that Lu = for almost every uj 

(u - P 2>go ,M« , u - P2, qo ,^) 2 qo41 > e 4Kr °C (u - P-i^u , u - P-i )qo ^u) 2 qQ pL ^ 

> C{u- P_ llWJ)M u , u - P-i m ,n u ) -i, qo ,v ■ 

Lemma 3.8. For every K > 1 there exists c = c(K) > such that for u G L 2 such that 
L u = for almost every uj and 

( u > u )-i, qo ,„ > c{P 2>qo> ^u, P 2>qo ^u) 2qo41 . (3.39) 

Proof of Proposition \2.1\ Of course Proposition 13.41 and Lemma 13.71 imply directly the 
spectral gap inequality for the Dirichlet form: 

S(u) ^ c K C (u - P_ 1)g0)M ii , u - P-i m ^u) _ 1)90)M for all it G J?~^ . (3.40) 

We now prove the self-adjoint property of A. It is sufficient to prove that the range of 
1 — A is dense in H' 1 (see [H p. 113]). For u,v G -D(/l), we have 

/R./S \Jo 9o/ ^JrJsQo 

v(6,u)J*u(6,u/)d0dn® n. (3.41) 



/a/(s) 2 

The right side of this expression is still defined for «,!)£ L 2 (A ® /u) (recall that A denotes 
the Lebesgue measure on S, and that we denote the usual scalar product on L 2 (A ® fi) by 
|| • \\2,fj,) and there exists c > such that 

(v, (l-A)u)_ lgo fi < c||u||2, M ||u||2,„, (3-42) 

Furthermore from (|3.40p and Lemma 13.81 we have 

(u, (1-A)u)_ ltg0jll > i||«||i t/x . (3.43) 

So the bilinear form (u, i-)- , (1 — A)u)_ l qQ is continuous and coercive on H' 1 x H~ x . 
If / G fl^" , the linear form v \— > (v , 9o „ is continuous on L 2 (A (X 1 /J,), therefore from 
Lax-Milgram Theorem we get that there exists a unique u G L 2 (A <%> /a) such that for all 
v G L 2 (A(g) /i) 

(«,(l-A)u)_ 1)ftJ>/1 = (v,f). lt90tlt . (3.44) 

Since 

<"' /> — =-/ R / 8 " ( ^ ) (r3 dM "' (s - 45) 

from (13.41 H we obtain that for almost and uj 

'«»>* + ^_/fe n ) (Si „ )d , = _r3>) d( .. (3.40) 



/„ ,0(9') 2<,„(9) ;„ V. ; > ' " J„ ,„(»') 

So it is clear that if / is continuous with respect to 6, then u has a version C 2 with respect 
to 9. Thus u G D(A) and applying de{qo{0)de-) to the both sides of this last expression, 
we get (1 — A)u = f . Since this kind of functions / is dense in H~ l , we can conclude that 
the range of 1 — A is dense, and that A is essentially self-adjoint. This completes the proof 
of Proposition 12.11 □ 
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4. Perturbation arguments (completion of the proof of Theorem 12.21) 

In this section we complete the proof of Theorem [221 Essentially, this section is devoted 
to computing the expansion of the speed c^(5) in. We first recall a lemma that gives a 
useful parametrization in the neighborhood of Mq. The proof of this lemma is given in 
[T5] , and it is used in the proof of Theorem 13.11 (see [TB] ) . 

Lemma 4.1. There exists a a > such that for all p in the neighborhood 

N a := UqeMoB^^^iq,^ , (4.1) 

of Mq there is one and only one q = v(p) £ Mq such that (p — q , doq) _ x qo = 0. Fur- 
thermore the mapping p i— > v(p) is in C 0O (Xj 1 ,Xj 1 ), and 

Dv(p) = P° v[p) . (4.2) 

Proof of Theorem \2.SX The existence and stability of a rotating solution q${6 — ip — Cn(S)t) 
of (jl.lOp (ip is arbitrary) has been established in Section [3] for 5 ^ do, see Theorem 13.11 
and the two remarks that follow it. We are left with proving Theorem 12. 2f 2). 

Thanks to the invariance by rotation, we can define qs such that v(qs) = qo. Now if we 
denote 

n s ■= q s -v (q s ) , (4.3) 
then ns verifies ns = 4>s(qo) and (see Lemma l4~T1) 



<A»*.fl&>_ w = 0. (4-5) 
Moreover the estimates we have on the mapping <fi$ in Theorem 13.11 give 

IK|| 2 ,„ ^ C5, (4.6) 

\\d e nsh,„ <CS. (4.7) 
Taking the derivative with respect to t, at time t = 0, we get (we recall the notation 

P ^\e,u) = qs(o-^-^(d)t)) ■ 

-c^5)(q' + d e ns) = d t pf . (4.8) 

So (jl.lOp at time t = becomes (recall that go is a stationary solution of (|1.12p ) : 

- c v( s ) (lo + dens) = An s - d e [n s {J * n 5 ) M ] - 5ujq' - 5ud e n 5 . (4.9) 
From (|4.6p we deduce the bound 

II^MJ*^]^ ^ \\Jhc 2 s 2 , (4.io) 

so by taking the scalar product of q' in (14. 9j) . using Q4.5j> . (j4.6|> . (14. 7p and the fact 

that L wdfj, = 0, we get that c^(5) is of order 5 2 . This implies, using the same arguments, 
that 

\\An s -6u>q \\_ ltfM = 0{5 2 ). (4.11) 

So 

\\A(n s - 5nW)\\^ = 0(5 2 ) , (4.12) 
and since ||(1 — Ap^'^uW-i^ ~ ||«||2,j* ( see <!3.42[) and (13.43P ). we have in particular 

\\ns-6nW\\2,v = 0{5 2 ). (4.13) 
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It allows us to make a second order expansion for c,(5) : taking again the scalar 
product of q' in (|4.9p . using the same bounds as for the first order expansion and (|4.12p . 
we get : 

c,(6) = 5 2 ^ j-r-r g °'-W + 0(j3) (4 . 14) 

Indeed, from flO^, \\iod e (n 5 - 5n^)\\_ h „ \\d 9 {(n 6 - 5n^)(J * n«) M ] ||0„[n«(J* 
(71,5 — ^n^ 1 )))^]!!-!^ are of order 5 2 and ||<96»[(^<5 — <5nW)(J* (n,s — ^n^ 1 )))^] ||_i iA1 of order 
<5 4 . Since c M (<5) is odd with respect to 5, the second order term in (|4.14p is equal to 0. It 
is possible to get this fact directly : we remark that satisfies : 

L qo [ n (1) d/i = f An {l) d{x = (I udfj) q' = 0, (4.15) 
Jr Jr \Jr J 

n^dfi,^) = U'Kq'o) i =0. (4.16) 



So since L qo is bijective on the orthogonal of q' in H~, (see |3j), we have J^n^ dfi = 

and (J * n^)fi = 0. On the other hand, since the operator A conserves the parity with 
respect to 8, n^ 1 ' is odd with respect to 8 and thus 

udenU ,q' ) = ! ! — L - / ) d9dn = 0. (4.17) 

Now back to (|4.9p : since 0,(5) is of order 5 3 and using J§ d/U = 0, we get 

A (n s - 8n,W - 6 2 u)d e n m ) = 0(5 3 ) , (4.18) 

and thus 

\\n s - 6nP> - 5 2 nW \\ 2j , = 0(5 3 ) . (4.19) 
This allows us this time to do a third order expansion in (14. 9p : 

This procedure may be repeated recursively at any order: we do not go through the details 
again, but we do report the result below (Remark 14. 2p and we point out that the 0(5 4 ) 
(I4.20p turns out to be 0(5 5 ), in agreement with the fact that 0,(5) is odd in 5. □ 

Remark 4.2. As anticipated above, one can get arbitrarily many terms in the formal 
series 0,(5) = 5^i=i 2 C2i+\5 2t+1 and the remainder, when the series is stopped at i = n, 
is 0(5 2l+3 ). In fact, by arguing like above, we have 

<a e [n( 2 )(J * n®),] + d 9 [nW(J * n^},} + wd e n^,q' )^ a< 

c 5 = T-7-7T ^ , (4.21) 

where 

An^ = d e [nU(J*nW),]+wd e nW- ^^^M^ q'^ (4.22) 

Wo>9o/-i,9o 

and 

An( 4 ) = ^[n^(J*n^) M ] + a e [n (1) (J*n( 3 )) M ] + ^n (3) . (4.23) 
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Actually, by induction we obtain 

C2, + i = : , (4.24) 

Wo,<7o)-l,go 

and 

= d e [n^(J*n^)^+wd e n^ i - 1 \ (4.25) 

k+l=2i,k>0,l>0 

n (2i+l) = J- do[nW(J * n( fc )) M ] + wd e n^ - c 2i+1 q' • (4.26) 

k+l=2i+l,k>0,l>0 

Since this procedure yields also tv**' for arbitrary j, one can generalizes also fj4. 19f) and, 
hence, (gHS])- 

5. Active rotators 

In this section we deal with the equation f|l . 13j) and we do it in a rather informal way, 
because on one hand a formal statement would be very close to Theorem 12.21 and, on the 
the other hand, the large scale behavior of disordered active rotators is qualitatively and 
quantitatively close to the non disordered case, treated in [8], in a way that we explain 
below. 

First of all, from a technical viewpoint the main difference between fjl . 13[) and (jl.4p is 
that (|1.13p is (in general) not rotation invariant, so the manifold Mg = {q^ + 4>(q^)} we 
get after perturbation is not necessarily a circle. Unlike Theorem 12. 2\ the motion on Mg 
is not uniform, and we describe the behaviour on Mg by the phase derivate ip. We follow 
the same procedure as in the previous section : if p\ is a solution f)1.13j) belonging to Mg, 
we define (see Lemma l4.ip 

m = v(p 5 t ) , and n\ = p\ - v(p S t ) . (5.1) 

In this context, (|4.9p becomes 



- + d t n\ = A^n\ - d e [n 5 t (J * nf) M ] - SU^ - 5Ud e n 5 t , (5.2) 

where A^ is the rotation of the operator A 

A^u(9,u) := ^Au(9,lj) - d e (q (6 - ?P)(J * u}^6) + u(9,lj)J * q (6 - ^ . (5.3) 

Note that we can reformulate the second term of the left hand side in (j5.2j) : 

d t nf = (5.4) 

So, as in the previous section, using the estimates on the mapping eft given in Theorem 
13.11 we get the bounds 

Hhv^CS, ||^|| 2)At ^ Ctf|$| and \\d e [n s t (J * nf)^ ^ \\J\\ 2 C 2 5 , (5.5) 

and we deduce the first order expansion 

Since tp is odd in 5 and the expansion can be pushed further in 5, this 0(5 2 ) is in reality 
a 0(<5 3 ) and one can actually improve this result both in the direction of obtaining a 
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regularity estimate on the 0(5 2 ) rest in (|5.6|) (like in [8] Th. 2.3]) and of going to higher 
orders (like in Remark I4.2p . 

However the evolution for small 5 is dominated by the leading order and from (|5.6p we 
can directly read that, to first order, the effect of the disorder is rather simple: in fact 



where c is such that J § (q^ - c) = 0, that is 1/c = 2TT(I (2Kr )) 2 (recall (pJ]) - (l23]l : 
this computation is analogous to (|4.15p ). Since the integrand depends on ui only via 
U, this integration can be performed first and the system behaves to leading order in 5 
as the non-disordered model with active rotator dynamics led by the deterministic force 
L U(-, uj)h( dw). The rich phenomenology connected to these models is worked out in [H 
Sec. 3]. 

6. Symmetric case: stability of the stationary solutions 

6.1. On the non-trivial stationary solutions (proof of Lemma 12. 3ft . We start by 
observing that in the case with no disorder the strict concavity of the fixed-point function 
has been proven in |1T|. Lemma 4, p. 315], in the apparently different context of classical 
XY-spin model (for a detailed discussion on the link with these models see [3]). We are 
going to obtain the concavity of ^ for small 5 via a perturbation argument, by relying 
on the result in |17j . 

Since ^fg is a smooth perturbation of one expects that the strict concavity of ^>q 
will be preserved to ^fg for small 5 > 0, namely sup x (*S>$)" (x) < 0. Nevertheless, an 
easy calculation shows that ^'(O) = 0; in that sense one has to treat the concavity in a 
neighborhood of as a special case. 

In what follows, we suppose that the coupling strength K is bounded above and below by 
fixed constants K m i n and K maiX : 

< K min < K ^ K m&x < oo. (6.1) 

We first prove the statement on the concavity in a neighborhood of 0: there exist rjo > 0, 
5 > such that for all K G [Kmim K max \, for all \i such that Supp(^) C [—1,1], ^ is 
strictly concave on [0, 770] • 

Indeed, one easily shows (using that the function x 1— > ^^(x) is odd) that we have the 
following Taylor's expansion: 

(#£)"0z) = -6D 5 (fi)K 3 x + e(x) , (6.2) 

where e(x) = o(x) as x — > and where for fixed fj,, we write 

D s (fi) := [ h(Su)fi(du), (6.3) 
Jr 

where 

j z \ 1 8w 2 . t . 

h{uj) '= 2(TT^) " (ITW ( } 

Note that the o(x) only depends on K max (in particular it can be chosen independently of 
/i). A closer look at the function h shows that there exists 5 > such that for all ji with 
Supp(;u) C [-1,1], D & (n) > \. If we choose r/ > such that i sup \e(x)\ < |-fC min 



' ?0 < x<r) 



then (^fg) (x) < for all < x < r/o, which is the desired result. 
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We are now left with proving concavity away from 0: namely, we prove that for all 
77 > 0, all i^ max , there exists Sq > such that for all K ^ K max , for all < 5 < 5q, for any 
measure \x such that Supp(/i) C [—1,1], ^ is strictly concave on [r],2K max ]. 

Indeed, using the strict concavity of \&o proved in |17] . there exists a constant a > 
such that for all x E [rj, 2K mscK ] , ^q(x) < —a < 0. But then, it easy to see that 

sup sup sup tf '(x)| 5o ->°0. (6.5) 

0<<5<<5 H, Supp(M)C[-l,l] xe[0,2K max ] 

If one chooses 5q such that the latter quantity is smaller than or equal to %, the result 
follows. The proof of Lemma 12.31 is therefore complete. □ 

6.2. On the linear stability of non-trivial stationary solutions. We now prove 
Theorem 12.51 along with a number of explicit estimates. 

Remark 6.1. Note that, since the whole operator is no longer self-adjoint nor sym- 
metric, its spectrum need not be real. In that extent, one has to deal in this section with 
the complexified versions of the scalar products defined in Section El (12. 8h and in Section 
O (|3.12p . Thus, we will assume for the rest of this section that we work with complex 
versions of these scalar products. The results concerning the operator A are obviously still 
valid, since A is symmetric and real. 

We will also use the following standard notations: for an operator F, we will denote by 
p(F) the set of all complex numbers A for which A — F is invertible, and by R(X, F) := 
(A — F)~~ , A G p(F) the resolvent of F. The spectrum of F will be denoted as cr{F). 

Decomposition of L„. In what follows, K > 1 and tq = ^o(2Kro) > are fixed. 

In order to study the spectral properties of the operator for general distribution of 
disorder, we decompose in (|2.19p into the sum of the self-adjoint operator A defined 
in (|2.6p and a perturbation B which will be considered to be small w.r.t. A, namely: 

Bu(9,uj) := -d e (u(e,u:)(J * s(q)}^ + e(q)(9,oj,S){J *u}^9) + 6uu(0,u)) , (6.6) 

where 

e{q) := (0, w, 6) ^ q(6, 5u) - q o (0), (6.7) 

is the difference between the stationary solution with disorder and the one without disor- 
der. 

Proposition 6.2. The (extension of the) operator A is the infinitesimal generator of a 
strongly continuous semi-group of contractions T^it) on H~^. 

Moreover, for every < a < ^ this semigroup can be extended to an analytic semigroup 
Ta(z) defined on A a := {z £ C ; | arg(z)| < a}. 

We recall here the result we use concerning analytic extensions of strongly continuous 
semigroups. Its proof can be found in [161 Th 5.2, p.61]. 

Proposition 6.3. Let T(t) a uniformly bounded strongly continuous semigroup, whose in- 
finitesimal generator F is such that G p(F) and let a G (0, The following statements 
are equivalent: 

(1) T(t) can be extended to an analytic semigroup in the sector A a = {A G C ; | arg(A)| < 
and ||r(^)|| is uniformly bounded in every closed sub-sector A' a , a' < a, of A a , 
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(2) There exists M > such that 

p(F) d E = {A G C ; | arg(A)| < - + a} U {0}, (6.8) 

anc? 

||i?(A,F)|| < ^, AgS,A/0. (6.9) 

Proof of Proposition ^.^ The proof in Section [3l Theorem 12. II of the self-adjointness of A 
shows that A satisfies the hypothesis of Lumer-Phillips Theorem (see |16t Th 4.3, p. 14]): 
A is the infinitesimal generator of a Cq semi-group of contractions denoted by Ta(£). 

The rest of the proof is devoted to show the existence of an analytic extension of this 
semigroup in a proper sector. We follow here the lines of the proof of Th 5.2, p. 61-62, 
in |16j . but with explicit estimates on the resolvent, in order to quantify properly the 
appropriate size of the perturbation. 

Let us first replace the operator A by a small perturbation: for all e > 0, let A £ := A—e, 
so that belongs to p(A e ). The operator A e has the following properties: as A, it generates 
a strongly continuous semigroup of operators (which is T^ j£ (t) = Tj^(t)e~ £t ). 

Since A is self-adjoint, it is easy to see that 

1 

|2KA)| 

and since the spectrum of A is negative, for every A G C such that 3?(A) > 0, 



VAeC\R,||i2(A,,4 e )||_ w ^ j^^, (6.10) 



For any a G (0, §), let 

E 

Let us prove that for A G E 



Q 



AeC; |arg(A)| < | + a} . (6.12) 



Note that (16331 is clear from (|6TT0l and (|6TTT1) when A is such that K(A) ^ 0. 
Let us consider a > 0, r G R to be chosen appropriately later. 

Let us write the following Taylor expansion for R(\,A £ ) around a + ir (at least well 
defined in a neighborhood of a + ir since a > 0): 

oo 

R(X, A e ) = ^R(a + ir, A £ ) n+1 ((a + ir) - X) n . (6.14) 

n=0 

From now, we fix A G E a with 3?(A) < 0. This series -R(A, A e ) is well defined in A if 
one can choose a, r and k G (0, 1) such that || R(o + ir, A £ ) qo |A — (<r + ir)\ ^ k < 1. 
In particular, using (|6.10|) . it suffices to have |A — (a + ir)\ ^ fc|r| and since a > is 
arbitrary, it suffices to find k G (0, 1) and r with |A — ir| ^ /c|r| to obtain the convergence 
of (gngp . For this A G Eqj with Sft(A) < 0, let us define A' and r as in Figured Then, 
| A — zt| ^ | A' — ir\ = sin(a)|r| with sin (a) G (0, 1). So the series converges for A G E a and 
one has, using again (|6.10|) . 

< (l-sin(o))|r| ^ l-L{a) 'W ^ 
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Figure 2. The set E Q . 

The fact that TA, s (t) can be extended to an analytic semigroup Ta e (z) on the domain 
A a is a simple application of (|6.15p and Proposition 16. 3( with M := 1 _ si 1 n( - a ) ■ 

Let us then define Ta(z) := e £Z TA )£ ( z ), f° r z £ A a so that Ta is an analytic extension 
of Ta (an argument of analyticity shows that Ta does not depend one). □ 

Remark 6.4. Note that estimate (|6.13j) is also valid in the limit as e — > 0: for all 

a 6 (0, § ), A € £„, 

Spectral properties of = A + B. In this part, we show that if the perturbation £> is 
small enough with respect to A, one has the same spectral properties for = A + B as 
for A In this extent, we recall that \x is of compact support in [—1, 1], and the disorder is 
rescaled by 5 > 0. 

Proposition 6.5. 

The operator B is A-bounded, in the sense that there exist explicit constants ok,s and bx,5, 
depending on K and 5 such that for all u in the domain of (the closure of) A 

II Bu H-W ^ aK > s I' U H-1,?o,m + bK > s I' Au H-l,«o,A» ' ( 6 ' 17 ) 
Moreover, for fixed K > 1, ax,s = 0(5) and bx,s = 0(5), as 5 — )• 0. 

The latter proposition is based on the fact that the difference e(q)(6,ui,5) = q(9,5uj) — 
qo(@) in (|6.7I) is small if the scale parameter 5 tend to 0: 

Lemma 6.6. For 5 > 0, let us define 

||e(g)IL:= sup \e(q)(8,u;,u)\. (6.18) 

0eS,H < i 

0<u<5 

Then for all K > 1, || e(g) = 0(5), as 5 — > 0. More precisely, for K > 1, 5 > 0, the 
following inequality holds: 

IK?) Hoc < ejf,*, (6-19) 
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where the constant £k& can be chosen explicitly in terms of K and 5: 

£K S .- 5 _ e ^ ( 2 + 3e 4 ^) e 14 ^ (l + 2vre 2 ^) , (6.20) 



7T 

where we recall that fg = ma,x(rQ,r$). 

Proof of Lemma \6.6l Recall that the disordered stationary solution q (|1.5j) is given by 

S(9,5u,2Kr 5 ) 

q{9 ' 5u}) '= Z(6 U ,2Kr s ) ' (6 ' 21) 
where S(9,lj,x) is defined in fll.6|) and that the non-disordered one (I2.1|) is given by 

?o(0) = S z(o,2Kro) = j e^oc° s S (fl) d e - Since <?(M W ) = -<Mi it suffices to consider 

the case <5w > 0. A simple computation shows that 

Z(8oj,2Kr 5 ) > 47r 2 e- iKrs e-^ s , (6.22) 

and that 

|5(0,O)| < 2ire 4Kr ° . (6.23) 
Using \q(9,5u) - q (9)\ < (2(O)|5(0,<M - 5(0, 0)| + |5(0, 0)||Z(0) - 2(<5o;)|), 

one has to deal with, successively: 

• for fixed 9 G S, \S(9,Su) - 5(0, 0)| < 5 • sup H ^ x |-^5(0, <5o;)|. A long calcu- 
lation shows that the latter expression \-4-S(9,5u))\ can be bounded above by 
8ir 2 e 4Krs e 47r<5 (2 + 3e 47r<5 ) , that is, 

|5(0, 5u) - S{9, 0)| s£ 58Tr 2 e 4Krs e 4wS (2 + 3e 4m5 ) . (6.24) 

• Using \Z(Su) - Z(0)\ = \f s {S{9,5u) - 5(0,O))d0| and QSMD , one has directly: 

|Z(£w) - Z(0)| < 516vr 3 e 4 ^e 4m5 (2 + 3e 4m5 ) . (6.25) 
Putting together (JS22]), (jE25D , (jSMD and (IBT251) . one obtains the result. □ 
We are now in position to prove the A-boundedness of B: 

Proof of Proposition \6.5[ B is ^4-bounded: let us fix a u in the domain of the closure of 
A. Then we have || Bu ||_ 1 = || Bu || 2 go , where Bu is the appropriate primitive of Bu, 
namely: 

Bu(0,w) := - (u(9,oj)(J*e(q))^ + e(q){9,uj,5)(J*u) fl (9) +5uju(9,oj)) 

+ f f 1 V V / n (^K J * £ (^))m + e(g)(e,Lj,5){J * u)^9) + ^(0,o;) d „ 

\JsQoJ \Js qo(9) 

(6.26) 

One can easily shows that there exists a constant depending only on K > 1 and 
(5 > such that: 

IWI-W < (6-27) 
Indeed, an easy calculation shows that |(J * ^ 4ff || e(q) and that 

\(J*v)p(-)\ < iff / sin(- -(^) 2 g (^d^) ||m|| 2j?0i/1 



(6.S 

< if I / (RjCv) ) || u || 2>W)>/i = K\\u || 2>W))/1 . 
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So we have for all 0,u (recall that Zq is the normalization constant in (|2.1|) ): 

\Bu(0,u>)\ < (^\\e(q)\\ x + 6\u\)\u\ + 2K\\e(q)\\ oo \\u\\ 2tq0tlt 

|2\| (6.29) 



+ Z 1 (4K\\e( q )\\ oo + 8\u;\)( f l -?f) 



Hence, inequality (|6.27p is true for the following choice of c$ s (recall that £k,s is defined 
in ([S^D ): 

i 

4] 5 := (6 (4Ke K! s + 5) 2 + 12K 2 Z 2 e^ 5 ) 5 . (6.30) 

Remark 6.7. Note that, thanks to Lemma 16.61 one has that = 0(5) as 5 — > 0. 
In order to complete the proof of the inequality (16.17p . it suffices to prove that there 

(2) (3) 

exist constants c K and c K , only depending on K such that, for all u: 

II U h qo ,» < c ? II ^« ll-l.go,/, + C K II u ll-W • ( 6 - 31 ) 

(2) (3) 

The rest of this first of the proof is devoted to find explicit expressions of c K and c K , 
and is based on an interpolation argument. 

For all integer n > 1, one can compute the linear operator / i— > f in terms of a sum of 
two integral operators, namely: 

/' = !„(/") + J n (/) , (6.32) 

where I n : f (-)■ J 27T i n (9 , ip) f (ip) dip (resp. J n : / i-> J Q 27r j n (9,ip)f(ip) dp) is the integral 
operator whose kernel i n (6,ip) (resp. j n (9,ip)) is defined by: 

' := j n (0,p) := - "("ty" 1 , < p < (9 < 27T, 

(6.33) 

k W¥>) := 2 V(2^ . Jn(0,¥>) := 2 7r(L-g)" ' < < V < 27T - 

Equality (16.32P can be easily verified by integrations by parts. Since, 

C\i n (e,v)\ dp < C\in(o,v)\ de ^ fa, 

(6.34) 

\ Jn (e, V )\ dip^^, C Me, v)\ de^ ^ , 

we see (cf. [121 p. 143-144]) that I n and J n are bounded operators on L 2 (E), namely: 

n — 1 7r(n — 1) 

So, applying relation (|6.32p for / = U we get, for /x-almost every u: 



This gives 



|W(#,w)| 2 d# . (6.36) 



2-7T ,, , n n(n + 1) „ , , „ ,„ 
u L„ «S u L + ^ 9„ ■ (6.37) 

"2,H ^ n _ 1 II H2,M 7t(tj - 1) " " 2 ^ V ' 
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Since II U L „ , , = II u \\ , „ „, it only remains to control II v! |L „ ,, with || Au || , , ,: like 
for the beginning of this proof for the operator B, we have || Au ||_ 1 = [| Au || 2 go „, 
where Au is the appropriate primitive of Au: 

Au(6,lu) := iu / (0, W )-( U (0, W )(J*( Z o)+goW(«/*w)^(0)) 

Using inequalities |(J*u)| M (-) < K^\\u\\ 2 ^, and / s J^li <; Z§e Kr ° (/ s |«(-, w) 2 |) *, an 
easy calculation shows that: 



|u'(-,w)| < 2|^n(-,w)| + 2Kr |u(-,o;)|+2 v ^^o(-)lhll2u + ^? e ^ ro ( / H'^f 



3? 



(6.39) 



and thus, 

||«'|| 2jM ^ 4Mn|| 2iAt + 4^(r 2 + 7rZ - 1 e 2 ^(l + 8r 2 ))^||n|| 2iM , (6.40) 

and by putting (|6.37p and (|6.4Q|) together we obtain 

8?r 2tt — * — i 



it 



n(n + 1) 
+ Tr(n-l) " ti "- 1 'W-' i 



(6.41) 



Let us choose the integer n 



16vr.Fr (r 2 + 7rZ^ 1 e Mfr »(l + 8r 2 )) 2 + 1 



so that 



-^4K (r 2 + vrVe 2 ^! + 8r 2 ))^ < I . (6.42) 
In this case, we obtain: 

e 2Kr ° 

4K" (r 2 + Tr^e 2 ^! + 8r 2 )) 2 

e 2^r Ag K ( r 2 + ^-^2^0 (! + g r 2^ I + 3^ 

+ " 1 — II « II > (6-43) 

16vr 2 if (r 2 + TrZ^etKroQ. + 8rg)) 2 

which is precisely the inequality (I6.3ip we wanted to prove. Inequalities (I6.27|) and (I6.31f) 

give the result, for clk,s '■= Crs " and bx,s '■= c£\ ■ . □ 

Proposition 6.8. For all K > 1, i/iere exists 5s(K) > suc/i that for all < 5 ^ ^(iC), 
i/ie operator is closable. In that case, its closure has the same domain as the closure 
of A. 

Proof. Let us choose £3 (if) > so that 

h KMK) < 1 ( 6 - 44 ) 
where 6^5 is the constant introduced in (|6.17p . then, for all < 5 ^ 5^(K), the operator 
B is ^-bounded with yl-bound strictly lower than 1. The result is then a consequence of 
Th. IV-1.1, p. 190 in [H]. □ 
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The spectrum of Lg . We divide our study into two parts: the determination of the position 
of the spectrum within a sector and its position near 0. 



Position of the spectrum away from 0. We prove mainly that the perturbed operator L°i 
still generates an analytic semigroup of operators on an appropriate sector. An immediate 
corollary is the fact that the spectrum lies in a cone whose vertex is zero. 

We know (Proposition 16. 2j ) that for all < a < 5, A generates an analytic semigroup 
of operators on A a := {A G C ; | arg(A)| < a}. 

Proposition 6.9. For all K > 1, < a < ^ and e > 0, there exists 64 > (depending 
on a, K and e) such that for all < 5 < 64, the spectrum of = A + B lies within 
&e,a ■= {A G C ; § + a ^ arg(A) < ^ - a} U {A G C ; |A| < e}. Moreover, there exists 
a' G (0, ^) such that the operator still generates an analytic semigroup on A a r . 

Proof of Proposition RTPl Let < a < ^ be fixed. Following (|6.17p and using (|6.16p . one 
can easily deduce an estimate on the bounded, operator BR^X^JV), for A £ ^>a>- 

II BR(X, A)u || _ w ^ a K>5 || R(X, A)u \\_ Xm ,i* + ^ H Ai2 ( A ' H - W 

1 

^ a ^(l-sin(a))|A| 11 U (6.45) 

+ (l + q ; r ) II « ILi 0n „ • 

' \ 1 - sm(a) / ' qo '^ 

Let us fix e > and choose 5 so that: 

max ( 1 —- + l) , ) < 1 . (6.46) 

\ \1 - sm(a) / (1 — sin(a)j£ / 

Then, for A G S a such that |A| > e ^ i^sinfa) ' we have 

\\BR(X,A)u\\_ 1>qo ^ < 5 || u ||_ 1|9WI . (6.47) 



In particular, 1 - BR(\, A) is invertible with (1 - BR(X,A)) 1 < 2. A direct 

-1,90,^ 

calculation shows that 

(A - (A + S))" 1 = fl(A, A) (1 - Bi2(A, ^)) _1 . (6.48) 
One deduces the following estimates on the resolvent: for A G S a , |A| > e, 

||E(A,L£)|| ^ 2 . (6.49) 

11 y 11 ±,<?o,m (l — sm(aJ)|A| 

Estimate (I6.49P has two consequences: firstly, one deduces immediately that the spectrum 
cr(Lq) of Lq is contained in O e>a : 

<r{I%) c{agC;|+«^ arg(A) < ^ - oX U {A G C ; | A| < e} . (6.50) 

Secondly, (|6.49|) entails that generates an analytic semigroup of operators on an ap- 
propriate sector. Indeed, if one denotes by L^ £ := — e, one deduces from (I6.50P that 
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G p(L^ 2£ ) and that for all A G C with 3i(A) > (in particular, |A| < |A + 2e\) 
\\R(X,Ll 2e )\l 1>q ^ = \\R(\ + 2e,L" q )\\^ qo4i < ^— ^ , 

^ (l-sin(a))|A|- (6 - 51) 

Hence, using the same arguments of Taylor expansion as in the proof of Proposition l6,2l and 
applying Proposition I6.3J, one easily sees that 2e generates an analytic semigroup in a 

(a priori) smaller sector A a /, where a' G (0, |) can be chosen as a' := \ arctan ^ 1 S1 ^ n ( a ) ^ _ 

But if £^25 generates an analytic semigroup, so does Lg. □ 

Position of the spectrum near 0. Let us apply Proposition 16.91 for fixed K > 1, a G (0, ^), 
p G (0, 1) and e := p\x, where we recall that \k is the spectral gap between the eigenvalue 
for the non perturbed operator A and the rest of the spectrum a (A) \ {0}. Let &f a ■= 
{A G & £ ,a ; 9i(A) ^ 0} be the subset of £jQ which lies in the positive part of the complex 
plane (see Fig. [3]). In order to show the linear stability, one has to make sure that one 
can choose a perturbation B small enough so that no eigenvalue of A + B remains in the 
small set Qf a - 




Figure 3. The set 6 £ , a . 

Since Xk > 0, one can separate from the rest of the spectrum of A by a circle 
centered in with radius (^-)\k- The appropriate choice of e ensures that the interior 
of the disk delimited by contains Qf a (see Figure [3]) . 

The main argument is the following: by construction of is the only eigenvalue (with 
multiplicity 1) of the non-perturbed operator A lying in the interior of % *. A principle of 
local continuity of eigenvalues shows that, while adding a sufficiently small perturbation 
B to A, the interior of still contains exactly one eigenvalue (which is a priori close but 
not equal to 0) with the same multiplicity. 

But we already know that for the perturbed operator = A + B, is always an 
eigenvalue (since L^q' = 0). One can therefore conclude that, by uniqueness, is the 
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only element of the spectrum of within and is an eigenvalue with multiplicity 1. In 
particular, there is no element of the spectrum in the positive part of the complex plane. 

In order to quantify the appropriate size of the perturbation B, one has to have explicit 
estimates on the resolvent R(X, A) on the circle ^ . 

Lemma 6.10. There exists some explicit constant c%> = c<g(K,p) such that for all A G ^ , 

11^.^)11-1,^ < <*> ( 6 - 52 ) 
|| AR(X, A) \\_ lmjtl < 1 + ( \k-cv. (6.53) 



One can choose &g as max ^^py, := 



R((,A) ||_ lgoM < sup ^-^ < (6.55) 



Proof of Lemma \6.1(A Applying the spectral theorem (see |23|, Th. 3, p. 1192]) to the 
essentially self-adjoint operator A, there exists a spectral measure E vanishing on the 
complementary of the spectrum of A such that A = j^XdE(X). In that extent, one has 
for any £ G 

KK,A)=[*m. (0.54) 

In particular, for £ G 

\ea(A) l A _ CI ^ X K 

The estimation (16.53P is straightforward. □ 

We are now in position to apply our argument of local continuity of eigenvalues: Fol- 
lowing [121 Th III-6.17, p. 178], there exists a decomposition of the operator A according 
to H~^ = Hq@H' (in the sense that AH C H , AH' C H' and PV(A) CD (A), where 
P is the projection on Hq along H') in such a way that A restricted to Hq has spectrum 
{0} and A restricted to H' has spectrum a (A) \ {0}. 

Let us note that the dimension of Hq is 1, since the characteristic space of A in the 
eigenvalue is reduced to its kernel which is of dimension 1. 

Then, applying [TJl Th. IV-3.18, p. 214], and using Proposition 16. 5| we find that if one 
chooses 5 > 0, such that 

sup (a Kj5 || R(X, A) || _ l qQ>lx + b K , s || AR(X, A) \\_ hqo ,») < ( 6 - 56 ) 

then the perturbed operator is likewise decomposed according to H~ „ = Hq ® H', in 

such a way that dim(ffo) = dim(.Ho) = L an d that the spectrum of L°i is again separated 
in two parts by ^€ . But we already know that the characteristic space of the perturbed 
operator according to the eigenvalue is, at least, of dimension 1 (since L^q' = 0). 

We can conclude, that for such an 5 > 0, is the only eigenvalue in and that 
dim(ffo) = 1. 

Applying Lemma 16.101 we see that condition (|6.56p is satisfied if we choose 5 > so 
that: 

a K ,scv + b K>s ( 1 + ( ) X K oA < 1. (6.57) 



2 

In particular, in that case, the spectrum of is contained in 

A G C ; - + a < arg(A) ^y-a C{zeC; R(z) ^ 0} . (6.58) 
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Finally, the following proposition sums-up the sufficient conditions on 5 for the conclusions 
of Theorem I2J3] to be satisfied: 



Proposition 6.11. Recall the definitions of ax,s a nd bx,s in Proposition \ 6.5[ If 5 > 
satisfies the following conditions 

or,8 < 1, 

4^,5 (- 1 - T ^ + l) < 1 , 

\ 1 — sm(a) / 

4a KtS (6.59) 



pXx (1 — sin(a)) 



a K /P- + b K ,S (l + l{p) ) < 1 ■ 



the conclusions of Theorem \2.5\ are true. 



Proof. One has simply to sum-up conditions (|6.44j) , (|6.46p with e = p\x and (|6.57p . (|6.60p 
can be obtained by (long) estimations on the coefficients a,K,5 and bx,5- □ 

Remark 6.12. The conditions in Proposition 16. Ill can be simplified. For example one 
can exhibit an explicit constant c such that if 5 satisfies 

,12** ^ m -20Kf f ™ QY (1 f l-rin(a) \ pA K (l - sin(a))e- 4g ^ 



5e iZ7TO < ce-" UA ^max 1, 



2 - sin(a) J ' K 2 

Xk 

K 2 e 4Kr si ^ + Xr ( i + fl+p) l{p) 



(6.60) 



the conditions in (|6.59p are fulfilled. Explicit estimates on the spectral gap \k can be 
found in [3l Sec. 2.5]. 
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Appendix A. Regularity in the non-linear Fokker-Planck equation 

The purpose of this section is to establish regularity properties of the solution of the 
non- linear equation ()1.13p (where we fix 5 = 1 for simplicity). Note that this case also 
captures the situation where U(-,uj) = oj (evolution (|1.4|) ). as well as the situation where 
U(-,-) = (evolution (|1.12[) ). In what follows we make the assumption that U is bounded 
and that for all u G Supp(^), 6 h-» U(6,uj) G C°°(S; R) with bounded derivatives. 

The existence and uniqueness in L 2 (A (8) ui) of a solution to (|1.13l) can be tackled using 
Banach fixed point arguments (see [IHl Section 4.7]), but one can obtain more regularity 
from the theory of fundamental solutions of parabolic equations. 
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More precisely, it is usual to interpret Equation (|1.13j) as the strong formulation of the 
weak equation (where v G C([0,T],A4i(S x R)) and F is any bounded function onSxI 
with twice bounded derivatives w.r.t. 9): 

[ F(O t u)v t (60, dw) = f F(9,uj)u {d9, duo) + ]- [ [ F"{0, u)v s (d0, dw) ds 

JRxS JRxS 2 ,/ Jr x § 

+ I I F'{9,uj)( I J(9--)dv s + U(6,uj))v s (de,du;)ds, 

JO JRxS \JrxS J 

(A.l) 

where the second marginal (w.r.t. to the disorder uj) of the initial condition z^o(d#, doj) is 
fi( duj) so that one can write 

z/ o (d0, dcj) = v$(de)n(doj), (A.2) 

where Vq is a probability measure on §, for (i-a.e. uj. 

As already mentioned, a proof of the existence of a solution on [0, T] of (jA.ip can be 
obtained from the almost-sure convergence of the empirical measure of the microscopic 
system |13| . One can also find a proof of uniqueness of such a solution relying on arguments 
introduced in |14| . 

The regularity result can be stated as follows: 

Proposition A.l. For all probability measure vo(d9, duj) = (d9)fi(duj) on § x 1R ; for 
all T > 0, there exists a unique solution v to (jA.ip in C([0,T],A4i(§ x R)) such that for 
all F G C(S x R), 

lim f F(e,u))v t (d9,du))= I F(9,u)i%(d9)n(du)). (A.3) 

*\0./]Rx§ JRxS 

Moreover, for all t > 0, vt is absolutely continuous with respect to Ai ® fi and for [i-a.e. 
uj G Supp(fi), its density (t, 9, uj) \— > Pt(9, uj) is strictly positive on (0, T] x S, is C°° in (t, 9) 
and solves the Fokker-Planck equation ()1.13[) . 

Proof of Proposition \A.1\ Let us fix T > 0, uj G Supp(^) and t i— >■ u t the unique solution 
inC([0,T],7Wi(SxR)) to (IA~T|) . Let us define R(t,0,u) := f RxS J(0--)di/ t + U(0,w) and 
consider the linear equation 

d t pt(9,io) = ^A Pt (9,uj)-d e (p t (9,uj)R(t,9,u;)}, (A.4) 

such that for fi-a.e. lj, for all F G C(S), 

/ F(9)pt(6,u>)d0 ^ / F(0)<(d0). (A.5) 

For fixed u; G Supp(^x), R(-, -,uj) is continuous in time and C°° in 9. 

Suppose for a moment that we have found a weak solution p t (9,u) to (|A.4p - (|A.5p such 
that for /i-a.e. uj, pt(-,uj) is strictly positive on (0,T] x S. In particular for such a solution 
p, the quantity L P t(9,uj) d9 is conserved for t > 0, so that p t {-,ui) is indeed a probability 
density for all t > 0. Then both probability measures ft(d9, duj) and p t {9,ui) d9fi(duj) 
solve 
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I F(9,uj)vt(d9, dco) = f F(9, oj)v ( d9, doj) + \ [ [ F"{6, u)v s ( d9, doj) ds 

JrxS JrxS * Jo JRxS 

+ [ f F'(9,u)R(t,9,uj)u s (d9, du)ds. (A.6) 

Jo JrxS 

By [H] or |144 Lemma 10], uniqueness in (jA.ip is precisely a consequence of uniqueness in 
(LATHI) . Hence, by uniqueness in (|A.6p . Vt{d9, du) = pt(9, w) d9/i( dco), which is the result. 
So it suffices to exhibit a weak solution p t (9,u) to ()A.4p such that (jA.5|) is satisfied. 

This fact can be deduced from standard results for uniform parabolic PDEs (see [2] 
and [6] for precise definitions). In particular, a usual result, which can be found in [21 §7 
p. 658], states that (|A.4p admits a fundamental solution T(9,t;9' , s,u>) (t > s), which is 
bounded above and below (see [21 Th.7, p.661]): 

. :ew ( -c(e-er \ <m ^ 8iU)< o =gxp r-^-^\ (A7) 



Note that the constant C > only depends on T and the structure of the linear operator 
in (|AH) (see % Th.7, p. 661] and [21 §1, p. 615]). In particular, since (6,u) h-> 11(9, lu) is 
bounded, this constant does not depend on u. 

Note that the proof given in [2j is done for S 6 1 but can be readily adapted to our 
case (0 e §). 

Moreover, thanks to Corollary 12.1, p. 690 in [2], the following expression of pt(9,u) 

Pt (9, uj)= f T(9, t; 9', 0, u)v$(d6') (A.8) 



defines a weak solution of (|A.4p on (0, T] x § (namely a weak solution on (r, T] x S, for 
all < r < T) such that (]A.5P is satisfied. The positivity and boundedness of p t (-,(v) for 
t > is an easy consequence of (|A.7[) . The smoothness of p.(-,uj) on (0, T] x S can be 
derived by standard bootstrap methods. □ 

We focus now on the regularity of the solution p t (9,co) of ()1.13|) with respect to the 
disorder u. We assume here that the initial condition vq is such that for all uj G SuppQn), 
Uq ( d9) is absolutely continuous with respect to the Lebesgue measure Ai on S: there exists 
a positive integrable function j(-,uj) of integral 1 on S such that Vq( dff) = 7(6*, uS) d9. Then 
we have 

Lemma A. 2 (Regularity w.r.t. the disorder). For every (io,#o) £ (0,oo) x S ; for every 
coq which is an accumulation point in Supp(fi) such that the following holds 

|7(0, oj) - 7(0, w ) I d# — > 0, as uj — > uj , (A. 9) 

then the solution p of (|1.13p defined on (0, 00) x§x Supp(/u) is continuous at the point 

(to,9o,LJo). 

Proof of Lemma \A.S\ . For any uj in the support of /i, let for all t > 0, 9 € 8 

u(t,9,oj) :=p t (9,u) -p t (9, uq), (A.10) 
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where (pt (•,•)) t > o is the unique solution of (|A.4|) . It is easy to see that u is a strong 
solution to the following PDE 



d t u(t,9,u) 



^Au(t, 6) - d e (u(t,9)R(t,6, coo)) 



n{t,e,u) 



(A.ll) 



where TZ(t,9,co) := dg\p t (9,uj) (R(t,9,co) — R(t, 0, coo))] an d with initial condition (since 
i#(d0) = -/(9,uj)d9 for all w) 

u(t J 0,a;)| tNl o = 7(0,a;)-7(0,a/o). (A.12) 

Then applying [61 Th. 12 p. 25], u(t,6,u) can be expressed as 

-t 



u(t,0,a;) = / r(0,t;0 / ,O,a;o)(7(^w)-7(0, W o))d0 / - 



o js 



For the first term of the RHS of ()A.13p . we have 



T(0, t; 9', s, uj )K{s, 9', u) d9' ds. 

(A.13) 



c_ 



| 7 (0,w)-7(0,wo)|d0' J (A.14) 



which converges to 0, for fixed t > 0, by hypothesis (|A.12p . 

Secondly, it is easy to see from the definition (|A.8|) of the density p and the estimates 
(|A.7p and [H Th.9 p. 263] concerning the fundamental solution T that both p t (9,uj) and 
dgpt(9,co) are bounded uniformly on (t,9,co) £ [0, T] x S x Supp(/i). In particular, a 
standard result shows that for fixed (t, 0), the second term of the RHS of (|A.13|) goes to 
as co — > ujo. But then the joint continuity of p at (to, w o) follows from (|A.8|) and uniform 
estimates on T (see [6j Th.9 p.263]). □ 
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